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The equilibrium conditions impose nontrivial geometrical constraints on the 
configurations that a two-dimensional foam can attain. In the first place, the three 
centers of the films that converge to a vertex have to be on a line, i.e. all vertices 
are aligned. Moreover an equilibrated foam must admit a reciprocal figure. This 
means that it must be possible to find a set of points Pi on the plane, one per 
bubble, such that the segments PiPj are normal to the corresponding foam films. 
It is furthermore shown that these constraints are equivalent to the requirement 
that the foam be a Sectional Multiplicative Vorono'f Partition (SMVP). A SMVP is 
a cut with a two-dimensional plane, of a three-dimensional Multiplicative Vorono'i 
Partition. Thus given an arbitrary equilibrated foam, we can always find point-like 
sources (one per bubble) in three dimensions that reproduce this foam as a gen- 
eralized Vorono'i partition. These sources are the only degrees of freedom that we 
need in oder to fully describe the foam. 



I. INTRODUCTION 

Cellular structures appear in a wide range of natural phenomena, and have puzzled and 

fascinated scientists for decades 0]. They can be generally described as packings of space-filling 
cells of roughly polygonal shape, separated by thin interfaces to which a surface energy is as- 
sociated. They arise as a consequence of competition between domains under the constraint of 
space-filling, and there is often some mechanism, such as the migration of a conserved quantity 
across the interface, which makes them evolve in time, i.e. coarsen. 

Foams such as those obtained by shaking soapy water are the examples of cellular structures clos- 
est to our everyday experience. Two-dimensional foams may be obtained by confining a soap froth 
between two closely spaced glass plates. In spite of their apparent simplicity, foams display much 
of the phenomenology appearing in coarsening cellular structures. Foams have been the subject of 
interest since their relevance in the problem of grain growth was pointed out by Smith || . Despite 
the attention they have received, the understanding of their dynamical properties proved to be a 
tricky problem. Even some of the most basic issues, such as their asymptotic scaling properties, 
has been a matter of debate until recently (see references above). Numerically exact simulation 
procedures as well as analytical J8j and numerical |^] approximations, have been useful in 
understanding the dynamics of ideal foams, but the system sizes accessible to available simulation 
procedures are strongly limited (for a complete account see refs. ||f|])- 

A satisfactory theoretical framework for the description of foams has not yet been achieved. The 
main result of this work consists in establishing a rigorous connection between foams and Vorono'i 
partitions(VP). This connection provides a set of fundamental degrees of freedom for the foam 
(the source's locations in space) and therefore constitutes a step towards the above mentioned 
theoretical understanding. On the other hand, an efficient method for the numerical simulation of 
ideal foams would also be highly desirable, and we propose that such a method could be obtained 
by exploiting the equivalence between foams and VP reported in this work. More precisely, we 
establish the following correspondence between equilibrated two-dimensional foams (EF) and a 
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generalization of VPs, the Sectional Multiplicative Voronoi Partition (SMVP): 

Given an arbitrary EF, it is always possible to find sources in space such that a SMVP with respect 

to them exactly gives this foam. 

In Section |l| the VP and its generalizations are reviewed. The simplest generalizations of the 
VP concept are the Sectional Voronoi Partition and the Multiplicative Voronoi Partition. These 
correspond to adding a constant to the (square) distance and to multiplying it by a constant 
respectively. At the end of this section a combination of both is introduced, the Sectional Multi- 
plicative Voronoi Partition. We will use this partition in order to describe foams. 
The demonstration of the above mentioned correspondence between foams and SMVPs is divided 
in two parts (Sections III and [V ) for clarity. In Section III , the recognition problem for SMVPs is 
solved. The recognition problem consists in giving the sufficient geometrical conditions that a cir- 
cular partition has to satisfy in order to be a SMVP. We will see that if a circular partition has all 
of its vertices aligned and admits an oriented reciprocal figure, then it is a SMVP. In other words, 
it is always possible to find sources in three-dimensional space that give this circular partition as 
a SMVP. In this section we also describe the procedure to find the sources when these conditions 
are met. Some material that is needed in this section is described in the appendix. 
Section ^ deals with the equilibrium conditions for foams and its geometrical consequences. We 
start by wri ting the f orce a nd pressure equilibrium conditions in compact form in Section 



IVA. 



In Sections pV B and IV C it is shown that an equilibrated foam has all of its vertices aligned 
and admits an oriented reciprocal figure respectively. Therefore all equilibrated foams satisfy the 
conditions required in Section III, and can be described as SMVPs. 

Finally in Section |y| the implications of this result are discussed, and some perspectives for future 
work are advanced. 



II. SPACE PARTITIONS 



A. Voronoi Partition 



Given a set of N point-like sources {fi} in n-space, the Voronoi Partition (VP) or tessellation 
of space with respect to {fi} is a classification of space into cells fii defined such that x £ if x 
is closer to fi than to any other source |hJ . 

n l = {xeR n / d(x,i) < d(x,j) Vj^i} (1) 

where d{x, i) denotes the distance between x and source fi. This construction is also known under 
the names: Wigner-Seitz cells, Dirichlet Tessellations, Thiessen Polygons. In two dimensions, the 
dual lattice of a Voronoi partition is called Delaunay triangulation. Cell fli can be seen as the 
"region of influence" of fi, were the sources competing for some spatially distributed resource. It 
is common to use these constructions, and their generalizations, as approximate models for cellular 
structures occurring in nature [^o|- |l2| , ^2| . 

Two neighboring cells fli and flj are delimited by an interface Tij of points x equidistant from fi 
and fj. 

T tj ={xe FT/ d(x, i) = d{x, j) < d{x, k) Vfc ^ i, j} (2) 

These interfaces are (n — l)-dimensional hyperplanes normal to fifj. In two dimensions f2i are 
convex polygons and the interfaces are straight lines (Fig. |l|). Three interfaces T^, Fjk and Tki 
meet at a vertex v^, which is equidistant from fi, fj and f^. This means that Vijk is the center of 
a circle through fi , fj and ff. . Vertices of higher multiplicity are possible for particular locations 
of the sources. For example a fourfold vertex would exist if four sources lay on the same circle. We 
will ignore the existence of these particular configurations, that is we will assume generic source 
locations. Under this assumption, all vertices are of multiplicity three. 

The concept defining Voronoi partitions is equidistance. A simple way to generalize them is 
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changing the way in which distances to the sources are measured. One lets each source measure 
the distance to points x according to its own "rule" . The interface is then the set of points 
x for which two neighbors "claim" the same distance. Two simple ways to do this are: adding 
an arbitrary constant to the square of the distance (Sectional Voronoi partition [ 10|, 13|JT9|] ) , and, 
multiplying the distance by a constant (Multiplicative Voronoi partition [jlO 141£^|2£|). 



B. Sectional Voronoi Partitions 

A Sectional Voronoi' Partition (SVP) is defined [ ^0|JT^Jl9[ | as a fc-dimensional cut of a higher 
dimensional Voronoi partition. The sources {Pi} of the sectional partition are defined as the 
projections of the original sources {fi} onto this lower dimensional fc-space. Cell VLi associated to 
source Pi is the intersection of cell Cli (associated to fi) with this fc-dimensional hyperplane. For 
example, take sources {fi} in three dimensions (Fig. ||) and construct a VP with three-dimensional 
cells Qi and plane interfaces Tij. Now project the sources fi = (xi,yi,Zi) onto the z = plane 
(which we call the Tl z plane) to obtain the projected sources Pi = (2^,1/j), which have associated 
"heights" Zj. Assign to each source Pi the intersection flj of f2, with IL Z . This defines the a 
sectional partition of IT with respect to sources {Pi} with heights Cells fii are thus defined as 

fii = {x e IT/ d 2 (x, i) + zf < (fix, j) + z) Vj ^ i} (3) 

where d{x, i) is the distance, on the H z plane, between Pi and point x. 
In the same way, interfaces are defined as 

r y ={a;elV d 2 {x,i) + z 2 = d 2 {x,j) + z 2 <d 2 {x,k) + z 2 Vk^iJ} (4) 

Properties of these partitions in two dimensions are: 

• Interfaces are straight lines and (generically) meet at triple vertices. ( In dimension d, 
vertices have multiplicity d+1). 

• Interface is normal to PiPj, but not in general equidistant from Pi and Pj. 

• The partition is unchanged if all Zi are changed according to: zf — -> z 2 + c 2 for c an arbitrary 
real. 

An important difference with the Voronoi partition, is that in this case there may be sources to 
which no cell is associated. This happens if (li is not cut by IT. It is easy to see that SVPs are 
equivalent to Laguerre Partitions ]l3[ ]. 

Appendix ^ discusses the recognition problem for this kind of partitions, that is giving the sufficient 
conditions for an arbitrary rectilinear partition to be a SVP. We will use similar concepts in 
Section |l| in order to solve the recognition problem for Sectional Multiplicative Voronoi Partitions. 
As an example of the application of sectional partitions, consider the case of sources fi whose 
location in space changes in time, giving rise to a time-dependent partition of II Z . If for example 
source f\ moves away from II Z , the two-dimensional cell Qi will shrink and finally disappear when 
fii no longer cuts II Z . Therefore the number of cells can vary without changing the number of 
sources. Thus sectional partitions can be used as dynamical models for crystal growth [^2| , and 
other processes in which some cells disappear or are created. 



C. Multiplicative Voronoi Partition 

The Multiplicative Voronoi Partition (MVP) JTJ,[l|,|26| is defined by assigning to each source 
fi a positive intensity aj, and defining the multiplicative distance d m {x,i) = d{x,i)/ai. The cell 
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fii associated to source fi is then the set of points x closer to fi (in terms of this multiplicative 
distance) than to any other source 

n l = {xeR n / ^£A < ^iA V jV<} (5) 

Interfaces are hypersphcrical surfaces (circle arcs in two dimensions) satisfying 

T lJ= {xeRV d ±A^ d ^lA < d ±2R vfc /i)i} (6) 

0>i CLj O-k 

In two dimensions, the circular interface of a MVP with two sources is one of the Apollonius circles 
(see for example pjfl) of those two points. In arbitrary dimensions let dy be the distance between 
two sources fi and /,-, and Ay = aj/ai- Without loss of generality we take Ay < 1. This means 
that fj has the smallest intensity and therefore £lj will be the interior of a hypersphere, while £li 
will be its exterior. The following properties are satisfied in any dimension. 

• fi is contained in f^. 

• Ty has radius Rij — Aydy/(1 — Ay), and its center Cy is on the straight line going through 
fi and fj. 

• Cy is located at a distance dy/(l — A?-) from that is, it never lays between /, and fj. 

Three sources give rise to a tessellation like the one shown in Fig. pi The exterior of the two 
"bubbles" is the infinite cell associated to the source with the largest intensity, f\ in this example. 
For some choices of {a^} the interfaces will not intersect. In this case the tessellation is simply 
a pair of circles, each containing one of the sources with smaller intensities, while the exterior of 
these two circles is the cell of the source with larger intensity. MVPs can be interpreted again as 
"areas of influence" of sources fi, but now each source has a different strength. 
Consider a MVP of two dimensional space with respect to a set of sources {fi}. If three interfaces 
Ty, Tjk and meet at a vertex (again, vertices of higher multiplicity are only possible for 
particular configurations, which we ignore here), then centers Cy, Cjk and C'ki lay on a line. 
The reason why this is so is simple. If two interfaces Ty and have a common point vyk 
, then their continuations must meet again at another point v*j k , the conjugate vertex. But 
then the third interface To- must also go through this point since d(v*,i)/a.i — d(v*,j)/a,j and 
d(v*,j)/dj = d{v*,k)/ak implies d(v*,i)/a,i = d(v*,k)/ak- Therefore two points exist (v and v*) 
at which all three circles intersect. Then the centers of these circles must be on a line (Fig. ||). 
Let us define a Circular Partition(CP) of two-dimensional space to be a classification of space in 
cells delimited by arbitrary circle arcs that meet at triple points called vertices. No three centers 
of these arbitrarily defined interfaces will in general be on a line. We will say that a vertex of a 
CP is aligned if the centers of the three interfaces defining it are on a line. As we saw, all vertices 
of a MVP are aligned in 2d. Therefore for each triple vertex of a MVP in 2d, the three sources 
and the three centers form the configuration ]23] (62, 4 3 ) of projective geometry. Since Cy cannot 
be between Pi and Pj , the three centers can only be on one of the two external segments of the 
configuration. 

In a three-dimensional MVP, four cells meet at each vertex Vyki , giving rise to six spherical 
interfaces (ij), (ik), (il), (jk), (jl), (kl). The centers of these interfaces are aligned in triplets so 
that the six centers also form the configuration (62, 43), this time in three-dimensional space. This 
has the implication that the six centers are necessarily on the same plane. 

The MVP was first introduced by Boots |L4| to describe areas of influence in geography. Ash 
and Bolker jl9| also studied the recognition problem, that is, under which conditions a CP is a 
MVP. The visual similarity between this kind of partitions and two-dimensional foams is striking, 
and suggests the idea to find a connection between them. Obviously the centers of the films of 
a foam must be aligned for each triple vertex, if the foam is to be described as a MVP, since 
MVPs are aligned. One finds |2(J that this alignment condition is indeed satisfied by all vertices 
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of arbitrary equilibrated foams. Despite this (which is a necessary but not sufficient condition 
for a foam to be a MVP), two-dimensional MVPs cannot describe all possible equilibrated foams 
in two dimensions. This work shows the reason of this limitation: In order to describe arbitrary 
equilibrated two-dimensional foams we must introduce the sectional variant of a multiplicative 
partition. In other words, instead of confining the sources to the plane, we let them exist in 
a three-dimensional space and obtain the foam as a two-dimensional cut of a three-dimensional 
MVP. 



D. Sectional Multiplicative Vorono'f Partition 

We will restrict the description to the case of a plane cut of a three dimensional MVP. The 
generalization to a fc-dimensional cut of an ra-dimensional MVP is straightforward. 
A plane cut, with a plane H x , of a three-dimensional MVP defines a Sectional Multiplicative 
Vorono'i Partition (SMVP) of U z (See Figs. | and |). The sources Pj of this SMVP arc the normal 
projections onto H z of the original sources /j, and may be seen to have as attributes both an 
intensity at and a height Zj. Cells Qj associated to sources Pi are defined as 

af aj 
In the same way, interfaces are circle arcs satisfying 

r<, = {* e n./ d ^±A = d Ih4±l < ^4±il vmm'} (8) 

", "j "f 

An example of a SMVP with three sources is shown in Figs. || and |[ We see that there are two 
vertices on II 2 at which the three interfaces meet. In a general case (for example in a partition with 
respect to many sources like in Fig. ^), if vertex exists, then the continuations of interfaces Ty, 
Tjk and Tki meet at a conjugate vertex v*j k also. This implies the alignment of centers Cij, Cjk 
and Cki, which could also be deduced from the alignment of E^, Ejk and E^i in three-dimensions 
(Fig. |). 

We notice that the SMVP is equivalent to a Multiplicative Laguerre Partition, since a SVP is 
equivalent to a Laguerre Partition [Q. As is usual in sectional partitions, in the SMVP there 
may be sources Pi with no associated cell, those for which the corresponding three-dimensional 
"bubble" f2j does not cut fl z . 



It was seen in section II C that the spherical interface Yij in three dimensions is cut normally by 



the segment fifj- As a projective consequence of this, the straight line containing Pi and Pj on 
II Z is normal to the circular interface Ty . In other words, Pj, Pj and are on the same line. 
Furthermore and since Eij never lays between /j and fj in 3d, we notice that Cij is always outside 
the segment P,Pj. 



III. WHEN A CIRCULAR PARTITION IS A SMVP 

In this section the sufficient conditions for a CP to be a SMVP are given. We first introduce 
the notion of oriented reciprocal figure for CPs, and then proceed to demonstrate that an aligned 
CP admitting such a reciprocal figure is a SMVP. 



A. Reciprocal figure of a circular partition 

We will now generalize the concept of reciprocal figure as appropriate for circular partitions. 
We say that a graph 1Z made of points Pi connected by edges (ij) forms a reciprocal figure for a 
CP V if, 
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points p are in correspondence with the cells tti of V . 
edges (ij) are in correspondence with the interfaces Tij of V . 

• for each in V , points Pi, Pj and lay on the same line. 

Consider two regions tti and f2j separated by a circular interface Tij. A reciprocal figure 71 will 
be said to respect orientation if for each Tij in V: 

• Cij is not between Pi and Pj. 



• starting from Cy and traveling along djPiPj to infinity points Pi and Pj are found in the 
same order as regions fli and Q,j are. 

We sa w alre ady that all SMVPs are aligned (all vertices satisfy the alignment condition introduced 



in Sec. II C). On the other hand, it is clear that the sources Pi of a SMVP form a reciprocal figure 
that satisfies orientation. Therefore 

V is a SMVP =>• V is aligned and admits an oriented reciprocal figure 

Figure Q shows four possible partitions that share the same centers C\i, C23 and C31. All four 
admit reciprocal figures. But only in cases a) and c) a reciprocal figure that satisfies orientation can 
be constructed. Therefore neither b) nor c) can be a SMVP. The reason is that the orientation 
condition cannot be satisfied if any of the vertices is not convex. A vertex is convex if all the 
internal angles formed by the interfaces are less than tt. Convexity of all vertices is clearly a 
necessary condition for a circular partition to be a SMVP. 

We are now ready to give the sufficient conditions for a circular partition to be a SMVP. We will 



show in Section III B that if a circular partition V is aligned and admits a reciprocal figure that 



satisfies orientation, then it is a SMVP. 

V is aligned and admits an oriented reciprocal figure V is a SMVP 



B. Finding the sources of a SMVP. 

Given a circular partition V of a region of two dimensional space, all whose (triple) vertices 
are aligned, we show here that, if V admits a reciprocal figure that satisfies orientation, then V 
is a SMVP with respect to sources fi located somewhere above points Pi. Our demonstration is 
constructive, that is we explicitly show how the sources and intensities are determined. We will 
use for this purpose the inversion transformation (Appendix |b|) to straighten each vertex in turn, 
that is, transform a circular vertex into a rectilinear one. This rectilinear vertex is one of a of 
a SVP. Sources are located in this straight representation (Appendix [a]) , and then transformed 
back to the original system. The intensities of the corresponding SMVP are obtained in this 
back-transformation, since the SMVP is invariant under inversion. 

Let v be an aligned vertex on which three interfaces {Fi2, T23, T32} meet, and {Pi, P2, P3} the 
three points of the reciprocal figure associated to the three bubbles sharing the vertex, as in Fig. ||. 
The conjugate vertex v* is obtained as the intersection point of the continuations of the interfaces, 
which happens at a unique point because of alignment. Let furthermore Xi be the straight line 
through v* and Pi. We start by showing that 

Theorem 1 A monoparametric family of triplets of circles {^12,^23,^31} exists, which has the 
following properties: 

• lo\2, W23 and W31 intercept each other on v* . 

• ^>ij -L T^ . 

• The intersection point qi between u>ij and Uik lays on \ . 
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The non-trivial content of the theorem is the fact that the three intersection points qi of these 
normal circles Uy lay on the lines Aj. For example take an arbitrary point q\ on Ai and draw a 
circle u>i2 through q\, v* and normal to T^. Let (72 be the intersection of u>i2 with A2. Draw now 
a circle W23 trough q2, v* and normal to 1^23. Let (73 be its intersection with A3. Then the circle 
W31 through v* , q\ and q% is normal to L13. 

Proof. The demonstration is done by first performing an inversion with center v* and arbitrary 
radius, whereupon the interfaces ry are transformed into straight lines Ty meeting at v' (Fig. |^). 
This transformed system will be referred to as the "straight" representation of that vertex. Primed 
names are used in the straight representation, with the exception of v*, whose original location is 
kept in Fig. [)] (v* itself is mapped to 00 by the inversion). Lines originally not going through v* 
are transformed into circles through v* , as is the case of {6*12, S23, Sm}, which form the reciprocal 
figure in the original system. Circles wy are transformed into straight lines wj - going through q[ 
and q'p respectively on Ai and Aj, which are invariant. 

The inversion transformation preserves angles, therefore circle <Sy is normal to the now straight 
interface Py. This means that its center Ly must lay on r' - . Consider the figure formed by v', 
L12, L23 and L31, and the six lines joining them. Maxwell J17| showed that a figure made of four 
points joined by six lines always has a reciprocal figure (to see this just consider the centers of the 
four circles going through three of these points). We will now identify it on Fig. I as the fi gurc 
formed by v*, q[, q' 2 and q' 3 . First we notice that circles SL and S'^ k intercept each other at two 
points v* and Pj, both on Xj. Therefore lines Ly-Zyfc are normal to Xj, and we have identified 
the first three lines of the reciprocal figure. The other three lines wy = q[q^ of the reciprocal 

figure have to be normal to the segments u£y , so they are normal to the interfaces Ty . Since the 
inversion preserves angles, this means that wy _L Ty in the original system (Fig. ||) . o 

We see that the existence of triplets of circles {u>i2 , ^Ti , ^31 } with the above mentioned properties 
is a consequence of the existence, in the straight representation, of a figure v*q' 1 q' 2 q' 3l which is the 
reciprocal of VL12L23L31. In other words, Theorem [l] means that in the straight representation, 
triangle q'iq 2 q' 3 forms a reciprocal figure for the rectilinear vertex. Clearly this reciprocal figure 
q'lq^q'z satisfies orientation if the original figure P1P2P3 satisfies orientation in the circular system. 
Therefore as discussed in Appendix |a|, it is always possible to find three sources f-, £ = 1,2,3, 
located at heights z[ above points q' i: such that this rectilinear vertex is a SVP with respect to 
them. Once these sources f[ are known in the straight representation, a second inversion around 
the same point v* takes us back to the original system and provides the original locations fi of 
the sources. 

If V123 is the first vertex to be considered in the system, we have one degree of freedom in the 
determination of the circles cjy, or equivalently of points q^. Once in the straight representation, 
there is one more degree of freedom: the determination of one of the heights z[. We will now show 
that, if this z[ is chosen appropriately, the three back-transformed sources fi are located above 
their respective points Pi. 

Consider a plane IL that is normal to II z and contains Ai, as shown in Fig. [l^. Draw a circle 
Pi through Qi, v* and normal to II Z . Let Ui be the intersection of this circle with the normal 
Q!j to IT Z through P^ Now under an inversion with center v* , point u[ is above q^, since circle 
(3i is transformed in a straight line /3^, normal to II Z . Straight line Ui normal to II Z through 
Pi is now a circle a[ through v*, and normal to II Z . Its intersection with II Z determines 
the image P[. Now we now that source f[ is located somewhere on (3[. The sources can be 
displaced vertically (simultaneously) according according to what we saw in Appendix [A], but not 
independently Fixing the height z[ of one of them determines the other two uniquely. What we 
want to demonstrate is that if one of the sources f[ coincides with its point u-, then the other two 
also do. In order to do this, take for example the sphere X\% through v*, /{, f 2 and normal to II Z . 
In the straight representation, interface T' 12 of the Vorono'i Partition with respect to sources /{ and 
f 2 is a plane, equidistant from /{ and f 2 , and normal to /{/^ This means that T' 12 is normal to 
A12. As a consequence, the intersection X\2(&, circle) of X12 with II Z is normal to the intersection 
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r' 12 (a straight line) of T' 12 with H z . Now assume that /{ = u[. This means that P{ £ X\ 2 since 
a[ is in this case contained in X\ 2 . But then the circle X\ 2 is the circle S[ 2 (Fig. ^) that goes 
through P[, P 2 an( i v * an d is normal to T' 12 . This implies that P 2 is also in X\ 2 , which in turn 
implies f 2 = u' 2 . The same reasoning can be repeated for the pair of sources 1 and 3. 
We have thus shown that all sources fi are, in the back-transformed, original system, located 
above the points Pi of the reciprocal figure. In other words, that the given reciprocal figure is the 
projection of the sources fi on the U z plane. 

Now we have to determine the heights Zi and intensities a^. The heights Zi of the sources can be 
found found by back-transformin g; th e sources f[ from the straight system. But in practice this is 
not necessary. By looking at Fig. |l(J we notice that triangle v* fiqi is rectangle. Therefore 

zf = \v*P\ x \P m \, (9) 

which suffices to locate the sources fi with the sole knowledge of points Pi and qi in the original 
system. Knowing now the spatial locations of the sources fi in three dimensions, the intensities 
a,i follow immediately since the vertex v* is equidistant [p9| from the three sources. This means 

a i =A\v*fi\ = A(\v*P i \x\v*q i \) 1 / 2 (10) 

Here A is an arbitrary positive constant. Alternatively we can get the same result by using the 
transformation properties (eq. |B^) of the dj's and the fact that the intensities are all equal in the 
straight system. 

We have used the inversion transformation to demonstrate that the given vertex is a SMVP with 
respect to sources fi located above Pi, and showed how the sources fi can be located without 
the need to explicitly perform an inversion for each vertex. All we need is to know the location 
of points qi, which we are free to choose for the first vertex under consideration. This first 



choice determines all subsequent qi points, and therefore all sources. It is clear from Fig. |10| that 
\v*qt\ > \v*Pi\ has to be satisfied. Therefore one has to choose the triplet of circles Wij such that 
this condition is verified for the three points qi. If one of the % coincides with Pi, this means that 
the corresponding height z, is zero, i.e. the source fi lays on the II 2 plane. A point qi closer to v* 
than the corresponding Pi is not acceptable. 

The general procedure to determine the sources can now be described. Assume we are given 
a circular tessellation which satisfies the required conditions of alignment and convexity, and 
that we are able to find, or are given, a reciprocal figure Pj for it. We start from one of the 
cells of the partition, as in Fig. O. Take an arbitrary vertex to start with, for example vertex 
an d det ermine tentative locations for point qo on Xq. This fixes q\ and q 2 as discussed in 



Section III B . After fixing qo the rest of the construction is uniquely determined. The locations 
<Z0i 9ij qi determine zq, z\ and z 2 , and the corresponding intensities do, a± and a 2 , through equations 
@ and (|h]). Now proceed to the neighboring vertex vq23- Two of the sources, /o and f 2 are already 
known, only the height 2:3 of f% has to be found. First one has to find the new locations of qo and 
q 2 when defined from vertex i>o 23 . I n order to do this, lines A° 23 from vertex tig 23 are drawn, and 
qo and q 2 are located using the fact that the heights z and z 2 are already known. Equation (||) 
implies 

\P r 3) \ = \PoC 2) \ l ^^\ (Ha) 
l>23)^ ol 

|P2 9 f 3) | = |P 29 f 2) |^l^ (Hb) 

Once the new qo and q 2 are known, the next step is to draw circles a->03 an d lj 2 ^ through them and 
normal to the respective interfaces. Their intersection gives the location of q% (this intersection 
always occurs on A3 23 , as Theorem |l| shows). This determines Z3 and a 3 using equations (||) and 
(|l0|). This procedure of triangulation is repeated until all sources are located. Eventually it may 
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happen that a point qi is found to be closer to the conjugate vertex v* than the corresponding Pj. 
This is not acceptable and means that the tentative starting position of qo has to be changed. It 
must be shifted away from Vq 12 and the whole procedure repeated. It is easy to see that taking a 
large enough value of q always solves this problem. 

One could ask whether this construction can be closed self-consistently. For example after de- 
termining /g from /o and in our example of Fig. one can go on with the procedure as if 
/i were unknown. Would the position of fx determined by /o and f$ be the same one as found 
initially? Alternatively: if we used fo and f± to determine /s instead of going around the bubble 
in the opposite sense, would its position be the same as found after going around the bubble? The 
answer is yes because of unicity. As discussed in Appendix [A|, the height of one of the sources 
of a vertex determines the other two uniquely. This means that /i,/2, • ■ ■ , /s are all uniquely 
determined by fo- 

IV. EQUILIBRATED FOAMS. 

In this section we show that a two dimensional foam in equilibrium satisfies all conditions 
required for it to be a SMVP. In order to do this let us first write down the equilibrium conditions 
for an arbitrary vertex of the foam, in compact form. We will consider the case of foams with 
arbitrary surface tensions, and also allow forces to act on the foam's vertices. 

A. Equilibrium equations 

Let v be the location of a vertex v at which three interfaces I\, T 2 and T3 meet, as shown in 
Fig. [l^. Each interface Ti is a circle arc with radius and center Cj. It produces on v, due to its 
surface tension , a force Tj of modulus 2r^ in the direction of the tangent to the film at v. These 
forces Ti can therefore be written as 

Ti = = (12) 

where K% — Ci — v and £j = 2ri/r,-. The asterisk stands here for a counterclockwise rotation in 
7r/2, so that e* = e y . 

Let us more generally assume that an external force F acts on v. Equilibrium of all forces acting 
on the vertex implies 

3 

F*+J2&Ki = Q t = 1,2,3 (13) 

i=l 

There is still one more condition, which is related to pressure equilibrium around the vertex. The 
pressure drop across an interface can be written as 

Ap 4 = 2^ = & (14) 
n 

We will adopt the convention for Ap to be positive if the pressure decreases when crossing the 
interface in the counterclockwise sense of rotation around v. This is of course related to a con- 
vention for the signs of the r^. In Fig. [l^ r\ and Api are negative according to this convention. 
Notice that pressure jumps and radii have different signs when considered from the two opposite 
vertices joined by a film. 

The fact that the total accumulated pressure drop around a vertex has to be zero implies then 

3 
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This is the pressure equilibrium condition for the vertex. Equations ( |13| ) and ( fig ) are satisfied if 
the vertex is equilibrated, and are together equivalent to 



3 

F*+J2^(Ci-x)=0, (16) 
i=l 

where x is an arbitrary point. This is what we will call the equilibrium condition for the vertex, 
and encloses both force and pressure equilibrium. 

B. Equilibrium implies alignement 

The alignment of the centers of a two-dimensional foam was already known long time ago by 
Plateau Eq] for equal surface tensions, and in the case of small self-standing clusters of 2 and 
3 bubbles. For the case of a cluster of two bubbles and zero load it is a trivial consequence of 
symmetry . The demonstration for the case of clusters of three bubbles is referred to by Boys 
as being "rather long and difficult" M|. 

It is not difficult to see that the alignment of the centers is in no way a property of clusters, and 
also not restricted to vertices with equal surface tensions and zero loads, but a general consequence 
of equilibrium. We will find that under very general conditions such as arbitrary surface tensions 
and external loads applied on the vertices, if a vertex is equilibrated then the centers of the three 
arcs converging to it lay on a line. 

Consider equation (|l6|), and assume for a moment that the external load is zero. Then C±,C2 
and C3 lay on the same line, as can be seen by taking x = C\ for example. Thus all vertices 
in equilibrium are aligned if no external force is applied on it. This alignment property is even 
true under non-zero load conditions, if the force is perpendicular to the line of centers. A load 
satisfying such condition will be called a normal load. The alignment condition has the geometrical 
consequence that the interfaces when continued, meet each other again at a unique point v* 
which we called the conjugate vertex [ p3| of v. This also means that the vertex v could be physically 
realized as a self-standing cluster of two bubbles (by continuing its interfaces) , and is a kind of a 
"separability" condition for the static equilibrium conditions, in the sense that each vertex of a 
foam might as well be that of an isolated cluster of two bubbles. Thus 

A vertex in equilibrium under a normal load is aligned. 

or ,equivalently 

A vertex in equilibrium under a normal load has a conjugate vertex. 

Next we would like to consider the existence of a reciprocal figure, since this is also a condition 
that has to be satisfied in order for a foam to be a SMVP. This condition must be separately 
considered. The reader may easily build examples of bubbles with n neighbors, all of whose n 
vertices are aligned, but yet do not admit a reciprocal figure. The reason is that the attempted 
reciprocal figure will not in general "close" around that bubble, the same case as described in 
Appendix 

C. Equilibrium implies Reciprocal Figure 

We will now show that if all vertices of a foam are equilibrated, then an oriented reciprocal 
figure exists for it. We start by considering a bubble and show that the reciprocal figure can be 
found for it. The figure for the whole foam can then be formed by patching together those of 
neighboring bubbles. Consider a bubble with n neighbors a = 1, 2, • • • , n, as shown in figure [l^. 
At each vertex v a of this bubble, three films T a , T a+ i and r" +1 meet. Interface r a separates 
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the central bubble from its neighbor a, while interface r" +1 is the limit between neighbors a and 
(a + 1). We will assume vertices v a to be in equilibrium under arbitrary normal loads F a . These 
loads we can generally write as 

F* = —R a +\C a +i + S a C a + T Q C" +1 (17) 

That is, we have decomposed the external load of vertex v a in the dependent basis formed by the 
three centers of the films meeting at the vertex. Because these centers are aligned, and the load 
is normal to the line of centers, the following condition is always satisfied by the coefficients. 

- R a+1 +S a +T a = 0. (18) 

The equilibrium condition ( |l6| ) now reads 

- + R a+ i)(C a+1 - x a ) + + S a )(C a - x a ) + (C +1 + T a )(C« +1 - x a ) = 0, a = 1, • • • n 

(19) 

and holds for arbitrary x a 's. The sign of £ Q is determined by the sign convention at vertex v a , 
therefore must appear with a minus sign in the equilibrium equation for vertex a. 
Since we use a dependent basis (equation. (|l7j)) the coefficients {R a +i, S a , T a } (the "representa- 
tion" of the load) are not uniquely determined for a given load. There is instead a monoparametric 
family of coefficients {R a +i, S a ,T a }, all satisfying both dl7| ) and (]l8|), for each load F a . We will 
use these degrees of freedom to choose a representation S a , T a } that satisfies 

S a = R a a = l,2,---,n (20) 

Notice that this condition relates the coefficients of two consecutive loads. We now show that such 
a representation always exists. We start by making the degree of freedom in the representation of 
the loads explicit. For arbitrary i/ a , we add the null vector ([l9|) multiplied by v a to the load (|l7j ) 
and get 

F* = —R a +iC a +i + S a C a + T a C^ +1 (21) 

where 

Ra+l = Ra+1 + V a (R a+ i + Ca+l) (22a) 

(R a +C a ) (22b) 

f a = T a + is a {T a + (22c) 

Condition (|2(]) then implies 

Rot ~ S a + i R a + ^a+l ,„„, 

Va + X = c T~c 7~c ' ( 23 ) 

which has always a solution in the generic case. 

Using this representation of the loads we can rewrite the equilibrium condition ( |l9| ) as 

-Q a +i(C a+1 -x a ) + Q a (C a -x a )+E a (C° +1 -x a ) = 0, a=l,...,n (24) 

where 

Q a = + R a (25a) 

E a = C +1 + T a (25b) 
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This amounts to hiding the loads in a redefinition of the surface tensions and pressures: £ Q — > Q a . 
As we have shown, this can always be done for aligned loads. We notice that the coefficients E a 
satisfy 

n 

^E a = 0, (26) 

Q = l 

as can be verified by subtracting Eq. (^J) with two different values of x a , and adding up the result 
for a — 1, ... ,n. 

More generally the fact that all vertices of the bubble are in equilibrium has the consequence that 

n 

Y / E a ((C^ 1 -x o ) = (27) 

a=l 

for Xq arbitrary. Condition ( pTj ) can be generalized to any closed path in the foam. The sum is in 
that case over all films cut by the closed path. The smallest such path, enclosing a vertex, gives 
the vertex equilibrium condition ( |l6| ) . Now that we have written the bubble equilibrium condition 
in the convenient form (p4|), we will write an algebraic condition for the existence of a reciprocal 
figure. 

The reciprocal figure was defined to be a set of points {Pq, Pi, ■ ■ ■ , P n }, each associated to a bubble 
(but not necessarily contained in it) , such that the straight line passing through Pi and Pj is normal 
to the interface between bubbles i and j. This means that Cij, Pi and Pj are on the same 
line, as in Fig. |l4|. We will require that Pq and P± be arbitrary (with the only condition that Pq, 
Pi are C\ are on a line) . This will allow us to patch together the reciprocal figures of neighboring 
bubbles to form that of the whole foam. This is equivalent to the translation and dilatation degrees 
of freedom existent in the definition of a reciprocal figure for a SVP. For arbitrary Pq, take Pi 
anywhere on the line (PqC\). All other points are now uniquely determined. P^ is located in the 
intersection of {PqC-i) and (P1C12), next P3 is found as the intersection of (P0C3) and (P2C23), 
etc. When the figure is closed with the last point P n , there is an extra condition, since it has to be 
the intersection of three lines: (PoC n ), (-P n -iC„_i„), (PiCi„). The construction of the reciprocal 
figure is thus an over determined problem, and would not in general be possible if the centers Cy 
were arbitrarily located. Let us now write the conditions for this reciprocal figure to close, in 
an algebraic form. The points Pi of the reciprocal figure are determined by a set of coefficients 
{A a ,B a } satisfying the following conditions. 

(P a - Pq) = A a (C a - P ) (28a) 
(P a+1 - P a ) = B a (C° +1 - P a ) (28b) 



Substituting 28a into f28q we get 



- A a+1 (C a+1 - P ) + Aa(l - B a )(C a - P ) + B a (C^ +L - P ) = (29) 



We will now show that the equilibrium conditions ( |24| ) ensure that (|2S|) always has a solution, and 
therefore that a reciprocal figure exists. Comparison with equation ( p4| ) lets us conclude that a 
solution will exist for Pq arbitrary, if we are able to find a set of coefficients {A a ,B a } that satisfy 



Q a+l Qa Pa 

A a +i A a (l — B a ) B a 

This is equivalent to 



H ai a = 1, • • • n (30) 



Qa+l — f-aAa+i (31a) 

Q a = fi a A a (l - B a ) (31b) 
E a = [L a B a (31c) 
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It is not difficult to see that these equations are satisfied if the /j, a are related by 

= Ha + E a (32) 
Starting from an arbitrary fj,i, this recursion relation gives us the following values of ji. Once all are 



known, equation (31a) provides the values of the A a , which in turn determine {P%, • • • , P n } using 
(28a). The condition that the figure can be closed is P a + n — Pa, and is equivalent to fj, a +n — fJ-a- 
This condition is satisfied because (|6|) holds, and therefore is a consequence of equilibrium. 
If we were given an arbitrary circular partition, it would not in general be possible to find a recip- 
rocal figure for it. The fact that this CP is an equilibrated foam imposes geometrical constraints 
on it ensuring, for example, that it admits a reciprocal figure. 

We have thus shown that, for an arbitrary equilibrated bubble, it is always possible to find a 
reciprocal figure. We can arbitrarily fix Pq since x a in ( p4| ) can be arbitrarily chosen, and we can 
also choose the "scale" |-Po-Pi| of the drawing at will, since the starting value /zi, that fixes this 
scale, is arbitrary. Therefore the reciprocal figures of neighboring bubbles can be patched together 
to form a reciprocal figure for the whole foam. Each selection of a starting point Pq and a "scale" 
\PqPi\ determines the other points uniquely, therefore there are three degrees of freedom in the 
determination of the reciprocal figure. Once Pq and Pi are chosen, all other points are found as 
intersections of two lines passing through the centers Cjj and one already existing point Pi . 



D. The orientation condition 



Now we have to show that it is always possible to find , am ong all possible reciprocal figures, at 
least one that satisfies orientation as defined in Section II D. In the first place, if all the surface 
tensions are positive then all films will be under traction and therefore the vertices will be convex. 
If there are non-convex vertices in the foam (which would happen if some of the films are com- 
pressed instead of stretched) then we know that it is not possible to satisfy orientation (Fig. 0). 
Positive surface tensions is thus a necessary condition |3(| for the foam to be a SMVP, although 
their modulus can be arbitrary for each film. 

The orientation condition could fail in the first place because a center lays in between two 
points Pi and Pj. It is always possible to avoid this by choosing Pi close enough to Pq. In this 
way all following points Pi are confined within a (arbitrarily chosen) small region of space where 
no center is located. This ensures that no center CV, lays between Pj an d Pj. 
Now regarding the second part of the orientation condition (Sec. II D), consider a vertex ^123, 
which is convex and aligned. Two orientations of the triangle P1P2P3 are possible, as shown in 
figures 0a and 0c. Notice that both constitute reciprocal figures for both vertices, but only one 
of them satisfies orientation in each case. In our construction of the reciprocal figure for the 
whole foam, we can decide the orientation of the initial triangle, choosing the one that respects 
orientation. The question is now if the correct orientation of this starting triangle ensures that of 
all subsequent ones, whose locations are determined by Pq and P\. To demonstrate that this is 
indeed the case, we notice that the triangle of figure 0a, if considered as a reciprocal figure for the 
vertex of 0c, has all three segments PiPj wrongly oriented. The point we want to make is that, if 
the vertex is convex, there are only two possibilities: either all pairs PiPj satisfy orientation, or 
all are wrongly oriented. Then if one of the pairs of a triangle forming part of a reciprocal figure is 
correctly oriented, the other two must necessarily also be. This demonstrates that if the starting 
pair PqPi is chosen with the correct orientation, then all subsequent triangles must be correctly 
oriented since they share at least a pair of sources with one preexisting triangle. Therefore in order 
to ensure correct orientation of the whole figure it is enough to correctly choose the orientation of 
the first pair. 
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V. DISCUSSION 



A dissection of space into cells separated by circular interfaces that meet at triple vertices is 
called a Circular Partition (CP). A two-dimensional foam therefore defines a circular partition of 
two-dimensional space. The equilibrium conditions for the foam impose geometrical constraints 
on this CP. We have here shown that the CP defined by an equilibrated foam is aligned and 
admits an oriented reciprocal figure. This result is valid in general for heterogeneous foams, each 
of whose films may have an arbitrary (positive) surface tension, and even if l oads a re applied on 



the vertices, with the sole requirement of equilibrium. We have seen in Sec. [II B, that any CP 



satisfying the conditions of alignement and existence of oriented reciprocal figure is a Sectional 
Multiplicative Voronoi Partition (SMVP). A SMVP is a plane cut of a Multiplicative partition, 
thus two-dimensional foams are plane cuts of three-dimensional foams, these last being a mul- 
tiplicative partition with respect to sources in three dimensions. Therefore given an arbitrary 
equilibrated two-dimensional foam, it is always possible to find sources {/$} in three-dimensional 
space, and amplitudes eij such that the given foam is a SMVP with respect to those sources. 
A first implication of this correspondence is the identification of a new set of degrees of freedom 
for the foam: the intensities ai and locations of the sources fi in three dimensions. This allows 
a more natural description of a foam, than the one that is done in terms of films and vertices. 
When a foam is interpreted as a tessellation of space with respect to some sources, we see the 
foam's films and vertices are secondary constructions, and their evolution is a consequence of that 
of the sources. The dynamical description is conceptually simpler using the SMVP interpretation. 
For example the processes of neighbor switching (Tl) and cell disappearance (T2) are described 
in a unified manner (Fig. |l5| ). Both are due to the fact that a fourfold vertex in 3d crosses the 
projection plane LT 2 . Depending on the spatial orientation of the vertex with respect to II Z , this 
is seen as a Tl or T2 process. 

An evolving foam can be seen as a particular instance of a Dynamical Random Lattice p7j, in 
which the evolution of a cellular structure is fixed by assigning a given dynamics to the sources 
of a mathematically defined tessellation. In the case of foams the dynamics is usually fixed by 
gas diffusion across the membranes. Alternatively other dynamical evolution rules may also be 
interesting, but in any case one has to make the translation to obtain the source's dynamics. The 
next step is then to find, for a given proposed dynamic evolution for the foam, the corresponding 
dynamics [^8| for the sources and intensities. 

Foams are usually studied inside a bounded region or cage, which imposes the only constraint that 
the films be normal to its boundary. Boundaries of this kind do not affect the fact that the foam 
is a SMVP. We have shown that the result holds with the sole requirement of vertex equilibrium. 
But it is obvious that SMVPs are always closed on themselves forming a self standing cluster, 
that is, there are no open film ends. This implies that even bounded foams must be a region 
of a larger self-standing cluster of bubbles that is closed on itself, and everywhere equilibrated. 
The point is not trivial in that it ensures that all films ending at the boundary can be continued, 
eventually forming new (phantom) vertices, and that the resulting foam will have all its vertices 
in equilibrium. We see then that there is no fundamental difference between bounded and self- 
standing foams, since all foams are regions of a self-standing cluster. This does not mean that the 
boundaries have no effect, which would of course be false. If the foam's dynamic inside the cage 
is for example produced by gas diffusion across the films, the evolution of the "phantom" bub- 
bles existing as continuations of the physical foam outside the cage, will not follow this dynamic, 
but a different one, which is determined by the constraint that the films be normal to the cage's 
boundaries. 

In the field of joint-bar structures, an old result due to Maxwell |l7|]l^] states that if a lattice 
accepts a reciprocal figure then it can support a self-stress, and conversely. More recently Ash 
and Bolker jl9| have shown that the existence of a reciprocal figure is sufficient condition for the 
lattice to be a Sectional Voronoi Partition. In this case there is the additional requirement that all 
vertices are convex, therefore all stresses in the lattice must be of the same sign, and the lattice can 
be an equilibrated spider web. A chain of results that span a century allow us to see equilibrated 
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spider webs as SVPs , and conversely. The alert reader may have noticed that equilibrated foams 
can be seen as a kind of "generalized" spider webs, in which the pressure difference between cells 
is the new ingredient, and is equilibrated by the curvature of the interfaces. It therefore turns 
out as no surprise that these generalized spider webs (foams) are equivalent to an appropriate 
generalization of SVPs, namely SMVPs, which include a multiplicative constant that gives rise to 
curved interfaces. 
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APPENDIX A: THE RECOGNITION PROBLEM FOR RECTILINEAR PARTITIONS 



A classification of space into cells VLi separated by rectilinear interfaces T^j is called a rectilinear 
partition of space. Given a rectilinear partition V of the plane, a reciprocal figure 7t(V)is a planar 
graph composed of sites i joined by edges (ij) satisfying |jl7|,[l9): 

• sites i of 1Z are in one-to-one correspondence with cells ili of V . 

• edges (ij) of 1Z are in one-to-one correspondence with interfaces of V. 

• edges (ij) of 7Z are normal to interfaces of V . 

Clearly a reciprocal figure is defined up to arbitrary global dilatations and translations, since 
angles are not changed by them. Therefore, if V admits a reciprocal figure, there are three degrees 
of freedom in its determination [ p~T|Jl~S| ] . An arbitrary partition V will not in general admit a 
reciprocal figure. We can see this with a simple example. Draw an arbitrary polygonal cell Slo 
with n faces, and take arbitrary rectilinear interfaces between its n neighbors (Fig. |l6| ). Now take 
an arbitrary point xq on the plane to start with, an assign it to the central cell. This starting 
point is arbitrary since a reciprocal figure is defined up to arbitrary translations. The other n 
points {xi, . . . , x n } associated to the external cells must be somewhere on the n rays stemming 
from xq and normal to the faces Toi of flo- The global length scale is also arbitrary so that say 
xqX\ can be freely chosen. Then we choose x\ somewhere on T\. Now point X2 is determined as 
the intersection of ray r 2 with the normal to face T 12 going through x\. This can be repeated to 
obtain all n external points, but in general the figure will not close, that is, segment x n x\ will not 
be normal to interface T\ ni whose orientation is arbitrary. 

The planar graph formed by joining the sources Pi of a SVP with edges (ij), one for each nonempty 
interface T^, constitutes a reciprocal figure for the SVP. Therefore every SVP has a reciprocal 
figure. 

V is a SVP => V admits a reciprocal figure 

Recently Ash and Bolker have shown that the existence of a reciprocal figure satisfying orientation 
fl9f is also a sufficient condition for V to be a SVP. A reciprocal figure TZ satisfies orientation if 
for each bond ij, the sites i and j are oriented in the same way as cells Qi and Qj are with respect 
to Tij. 

V admits an oriented reciprocal figure =>■ V is a SVP 

This result solves the recognition problem for SVPs. The orientation condition can only be satis- 
fied if all vertices of V are convex. We will say that a vertex is convex if the internal angles formed 
by the interfaces are all smaller than n. Figure |l^ shows two partitions with three cells. One of 
them is convex, the other is not. In the second case it is not possible to find a reciprocal figure 
that satisfies orientation, and therefore it cannot be a SVP. 

Given a rectilinear partition V* , and an oriented reciprocal figure for it (as in Fig. |l8| ), it is always 
possible [|l9| to find sources {fi} in three dimensions, located at heights z% above the points Pi, 
such that V* is the section with II Z of a three-dimensional VP with respect to {fi}. The procedure 
to determine the heights Zi can be easily described. Vertices Vijk of V* will be equidistant from 
sources fi, fj, fk (see Fig. [l9]). Start from an arbitrary vertex, say Vijk, and draw a spherical 
surface of arbitrary radius r^k = ro with center at that vertex. Now define sources fi, fj and 
fk as the intersections of this surface with the verticals (normals to IT Z ) through Pi, Pj and Pk 
respectively. Next go to a neighboring vertex, which shares two sources with this one. Let us call 
it Vij m . For this vertex, only source f m has to be located since fi and fj are known. Draw a 
spherical surface with center Vij m and containing fi and fj . Both will be simultaneously contained, 
since interface r^ , on which Vij m is located, is equidistant from fi and fj. The radius r,j- m of this 
surface is determined by the locations of fi and fj, which in turn is fixed by ro. Its intersection 
with the vertical n m through P m determines f m . If this spherical surface does not intersect n m , 
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just choose a larger value of ro and start all over again (from the initial vertex). The construction 
proceeds in this manner until all sources have been determined. As mentioned, the initial value 
of ro is tentative, in the sense that it may have to be modified (increased) if at some point during 
the procedure, a normal line is not cut by the corresponding spherical surface from the vertex. 
It is easy to see that increasing the value of the starting radius ro is always enough to solve this 
problem. 

There is thus one degree of freedom in this construction (ro). We conclude that, given a two- 
dimensional partition V that admits a reciprocal figure, there is a four-parametric family of source 
locations such that V is a SVP with respect to them. Three degrees of freedom come from the 
determination of the reciprocal figure itself (since a dilatation and/or translation of a reciprocal 
figure is again a reciprocal figure) and the last one from ro- This last degree of freedom results 
from the fact that a SVP is invariant if all heights are changed according to zf — ► zf + a 2 with a 
arbitrary (see equations @ and (Q)). 

Reciprocal figures were first studied by Maxwell 0,^8| in relation with the rigidity of bar-joint 
frameworks in the plane. Maxwell pointed out that frameworks that have a reciprocal figure are 
able to support a self-stress, and conversely. The reason is that the edges of the reciprocal figure 
can be taken to represent forces transmitted by the edges of the original framework (rotated by 
7r/2). Since these edges form closed polygons, the existence of a reciprocal figure implies the ex- 
istence of an equilibrated set of internal stresses in the absence of external load. The addition of 
the orientation condition (a condition not required by Maxwell's definition of a reciprocal figure) 
has the statical consequence that all signs of the stresses are equal, for example all traction or all 
compression. It is clear that no equilibrium is possible in the case of Fig. [17p, if all three stresses 
are to have the same sign. Fig. [T^ b on the other hand, can be in equilibrium under compression 
or traction on all three interfaces. 

The conclusion is that every SVP is an equilibrium configuration for a spider web Jl9| ] , and con- 
versely, each such equilibrium configuration is a SVP. 

The existence of a reciprocal figure has also projective consequences, which have been studied by 
Crapo fU and Whiteley glj. 



APPENDIX B: THE INVERSION TRANSFORMATION 

We briefly describe here a geometric transformation called inversion p4|] . We will find it ex- 
tremely useful for our purpose of discussing circular partitions. An inversion with radius k around 
a point O located at r*o transforms a point P located at r*o + f into a point P' at ro + r ' satisfying 

k 2 

r"=^ r (Bl) 

where r = \r\. The sphere of radius k and centered at O is invariant under this transformation, 
while the inside and outside of this sphere arc interchanged. Obviously this transformation is self- 
inverse. Let us now describe some of the properties of this transformation in two dimensions [ p4[ . 
Most of them apply trivially in higher dimensions. 

• Circles not through O are transformed in circles not trough O. 

• Circles through O are transformed into straight lines not through O. 

• Straight lines not through O are transformed in circles through O. 

• Straight lines through O are invariant. 

• Angles are preserved (in modulus) by the inversion. 
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Given two points Pi and Pi at distances r± and ri from the inversion center, the distance dn 
between them transforms as 



(B2) 



Using this result we can easily see that the inversion is a "symmetry" of a MVP in any dimension if 
the intensities are also appropriately transformed |L9| . More precisely, a MVP of R n with respect 
to sources {Pi} with intensities {a^} is transformed by an inversion into a MVP with respect to 
{P/} with intensities {a^}, where the new intensities satisfy 



a, 



A 



(B3) 



Here Aq is an arbitrary prefactor, the same for all Oj's, and r, is the distance between source i 
and the inversion center O. To see this it suffices to demonstrate that if x £ Ti,- then after an 



inversion, x' £ , which is easily done using (§), J32| ) and @. The inversion transformation is 
of course also a symmetry of the SMVP (Section IIP ), if the inversion center O is on the cutting 
plane IT Z , since in this case the inversion leaves this plane invariant. 

If the inversion center O happens to be located on an interface T^- of a MVP (initially a spherical 
surface), the transformed interface is a plane not through O. The resulting interface thus 
corresponds to a Voronoi Partition with respect to sources i and j in their new locations. Therefore 
the transformed intensities a! i and a' have to be equal after the inversion, which is verified using 

(H) 



0£r, 



(B4) 



In the same way a SMVP with respect to two sources i and j is transformed into a SVP if 



O £ Tij. The intensities ai are transformed according to (B3), where is the distance between O 



and so urc e ft in three-dimensional space. The way in which heights Zi transform is easily found 
using (Bl). Notice that if the inversion center coincides with the location of a conjugate vertex v* , 
then the transformed partition has a rectilinear vertex since three interfaces are simultaneously 
transformed into straight lines. We will use this prop erty of the inversion transformation when 
solving the recognition problem for SMVPs in Section [II B. 
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FIGURES 




FIG. 1. Shown is a Vorono'i partition of two-dimensional space with respect to point- like sources {/;} 
(black dots). Each source /, has an associated cell Q», which is the subset of space that is closer to /; than 
to any other source. Interfaces Yij (thick lines) are equidistant from the sources whose cells they delimit. 
Three interfaces F^-, Yjk and F^ meet at a vertex Vijk, which is therefore the center of a circle (dashed) 
trough the three corresponding sources. 




FIG. 2. A Sectional Vorono'i Partition (SVP) of the plane is defined by intersecting a three-dimensional 
Vorono'i Partition with that plane. Sources P; of this SVP (black dots) are the projections of the 
three-dimensional sources /; (white dots) onto the plane. Interfaces Yij between cells are rectilinear, 
normal to PiPj and meet at triple vertices Vij k , but are not equidistant from the sources. Sources with a 
too large value of z, (for example PiO in this figure) do not have an associated cell. 
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FIG. 3. A Multiplicative Vorono'i Partition with respect to three sources in two dimensions. The source 
with the largest intensity at (source /i in this figure) has an unbounded cell. The other two sources are 
within their associated cells. Angles formed by the interfaces at the vertices are depend on the positions of 
the sources and the values of the intensities. For certain values of {a,i} the interfaces may not intersect. In 
this case the partition would be composed of two disjoint circles containig one of the sources with smaller 
intensities each. The MVP with respect to three sources in three-dimensional space is a cluster of two 
spherical bubbles (Fig. ^, dashed lines). 
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FIG. 4. A Sectional Multiplicative Vorono'f Partition (SMVP) is defined as a plane section (black 
circles) of a three-dimensional Multiplicative Vorono'f Partition (dashed lines). Here an example for three 
sources. Sources Pi of the SMVP (on the II Z plane, white dots), are the projections of sources fi of the 
MVP in three dimensions (black dots). The centers Etj of the 3d MVP are aligned, which implies the 
alignement of the centers dj of the SMVP. Segments /»/,• are normal to spheres F;j, therefore PiPj are 
normal to circles on the plane. 
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FIG. 5. An upper view of the SMVP with respect to three sources defined in Fig. H. Notice that 
sources Pi are not all contained in their cell Qi, as happens in the case of the MVP (Fig. |[). 




FIG. 6. A SMVP of the plane with six sources. The projections Pi of the sources are indicated as 
small dots in the figure. For each one of the vertices in this partition, the three interfaces can be continued 
and they will meet each other again at a conjugatevertex. 
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FIG. 7. All these four circular vertices admit a reciprocal figure. Two of the vertices, a) and c), are 
convex and therefore a reciprocal figure satisfying orientation is possible. The other two are not convex 
so no reciprocal figure can satisfy orientation on them. 
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FIG. 8. A vertex v with its reciprocal figure P1P2P3, showing circles uiij. These circles are normal to 
the interfaces Tij, and intersect each other at points qi on Xi. 





FIG. 10. Shown is a plane normal to II z and containing A^. This figure shows how the height 2, of 
source /» = Ui in the original system is determined, given the positions of Pi and qi. Primed variables 
correspond to positions after the inversion, that is, in the "straight" system. Sources = u[ are above 
q'i, and therefore the back-transformed sources fi are above Pi. 
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FIG. 11. The points Pi of the reciprocal figure (empty circles) are the projections, on the xy plane, of 
the sources ft in three-dimensional space. The lines PiPj must thus be normal to interfaces Fij. Black 
dots indicate the location of the conjugate vertices v*j k for each vertex Vijk in the figure. 
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FIG. 12. A generic vertex v, at which three interfaces Fi, T2 and T3 meet. We regard pressure drops to 
be positive if the pressure decreases when the interface is crossed in the counterclockwise sense of rotation 
around v. Therefore in the case shown in this figure, and r-j, are positive, while n is negative. 




\ 

FIG. 13. A bubble surrounded by n neighbors a = 1, . . . , n, showing the naming convention for the 
interfaces. Notice that film curvatures are of opposite signs when considered from each of the two vertices 
at their ends. Thus for example £2 appears with + sign in the equilibrium condition of vertex V2, but with 
— sign in that of vi (See equation Il9[) . 
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FIG. 14. The same closed bubble of Fig. |l^, showing the film centers dj (empty dots) aligned in 
triplets, and the sites 0, 1, . . . , 6 (black dots) of the reciprocal figure. 




FIG. 15. A three-dimensional Multiplicative Vorono'i Partition gives rise to vertices which simultane- 
ously belong to four neighboring bubbles in a generic situation. Six spherical interfaces converge at these 
vertices. When one such vertex (white dot) crosses the projection plane II Z , depending on its orienta- 
tion this can result either in a) a T2 process (disappearance of a triangular bubble) or b) a Tl process 
(neighbor switching) for the two-dimensional foam. 
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FIG. 16. An arbitrary rectilinear partition will not in general admit a reciprocal figure. The one in 
this example does not admit a reciprocal figure, and therefore cannot be a Sectional Vorono'f Partition. 




FIG. 17. Both these partitions with three cells admit a reciprocal figure. The one on the left does 
not satisfy orientation and therefore cannot be part of a Sectional Vorono'f Partition. The one on the 
right satisfies orientation, and therefore it is a SVP. It is then possible to find sources /i, f% and fs in 
three-dimensional space, such that this partition is obtained as the cut z — of a three-dimensional 
Vorono'f Partition with respect to {/;}. Sources /, will be located at heights Zi above points {1,2,3} in 
this figure. 



29 



FIG. 18. A convex rectilinear partition (thick lines) with cells {A,B, ...,/} and its associated recip- 
rocal figure (thin lines) with sites {a, b, ■ ■ ■ ,i} (black dots). This reciprocal figure (or any other figure 
obtained from dilatation and translation of this one) satisfies orientation, and therefore the partition is a 
Sectional Vorono'i Partition. Sources of this SVP are located above points Pi. 




FIG. 19. Given an oriented reciprocal figure {Pi,P2,PS} (white dots) for a rectilinear partition 
(interfaces T^), the sources fi in three dimensional space can be located above the points Pi. These 
sources must be equidistant from the vertex, therefore they are at the intersections (black dots) of a 
spherical surface centered at the vertex «i23 with the normals through Pi. 
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